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AbstfCt 

U« pr«a«nt a th«nul Bodcl for th« lov«r which is coiustructod 

froB pstrologicslly dsrlvcd •■tlMtss of th« toeq^oratur* in the tcanaltion 
sona and fron an adiabat baaed on the theraal properties of NgO and Si02 
neaaured at high preaaurea. Superadlabatic contributions to the geothem 
through the lower mantle, including thoae from pr^:^aible phaae tranafomationa 
but excluding thoae from theraal boundary layera, are negligibly amall. 

A thermal boundary layer is required at the base of the mantle in order 
to aatiafy our estimate of the lowest possible temperature in the core 
(about 2800K); its thickness of about 100 km is well constrained. Even 
so, our model suggests that we either require a significantly larger heat 
flux from the core than has been considered reasonable (> 50 mWa"^) or 
we require the presence or a further thermal boundary layer within the 
lower mantle in order to arrive at more likely core temperatures (about 
3200*3500K). This latter alternative appears the more plausible, and 
such a boundary layer requires a barrier to convection. Guided by the 
present seismological evidence, we suggest that a thermal boundary layer 
is associated with a chemical discontinuity either at the top of the 
lower mantle or near its base (D" region). If seismological data require the 
D" layer to be significantly thicker than about 100 km, this is support 
for a chemically distinct D" region with multiple thermal boundary layers. 
Otherwise, the results of this study suggest that the upper and lower 
mantle nay be chemically diatinct, hence precluding whole-mantle convection. 
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Introduction 

Thu conposltlon and dynamic atata of tha aacth'a lovar mantla hava 
avokad conalii!:iifabla attention recantly> In particular with raferanca to 
arguaanta balng preaanted for or aglnat nuintla-wlda convection [Davlaa 
1977; Jordan* 1977; O'CoanalX* 1977; Richter and McKansla* 1978]. 

Slallarly* It haa bean argued that aalanologlcal data* coablnad with 
Masured and v^strapolatad physical properties of candidate minerals* do 
or do not require a difference In bulk chemical composition between 
tha upper and lower mantle [e.g. Birch* 1961; Press* 1968; Anderson* 1968* 
1977; Anderson* at al. 1972; Mao* 1974; Davies* 1974; Rlngwood* 1975; 
Burdick and Anderson* 1975; V/att, at al. 1975; Liu* 1979]. These two 
questions are* of course* coupled In that whole'-mantle convection Is not 
compatible with a compositional change between the upper and lower mantle. 

We propose a different approach for addressing this problem by 
examining the thermal state of the lower mantle, and the ensul.ng Impli^ 
cations for convection and chemical composition throughout the mantle. 

We find that the temperature gradient through the lover mantle Is well 
enough constrained to Indicate the existance of a thermal boundary layer at 
the base of the mantle. Furthermore* our preferred solution requires a 
second boundary layer elsewhere within the lower mantle* most likely 
either close to Its base or associated with the transition zone. This 
strongly suggests the existence of a chemical transition which Is a 
barrier to convection either at the top of the lower mantle or in a complex 
region at Its base. A thermal structure with only one thermal boundary 
layer within the lower mantle (presumably at the core-mantle boundary) 

Is only acceptable If the outer core Is at a lower temperature than has 
previously been estimated or if the heat flux from the core Is significantly 
larger than has been considered reasonable. 
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Adiabatic Gr«dl<nf 

Wa axaaina tha follovlnt aquation fo\r tha Incraaaa in taaparatura 
froK tha top of tha lowar Bantla (T(700)) to tha top of tha outar cora 
<T(2900)>: 

T(2900) •' T(700) - lATg + 0Azl + E AT(P.T.) + E AT(B.L.). lU 

Tha brackatad tan contalna tha naarly adiabatic taaparatura riaa through 
the lower aantlUf appropriate to a convacting and hoaogeneoua aystem 

(6 ia the average auperadiabatlc gradient and Az la the thickneaa of the 
lower aantle). Tha remaining tana allow for temperature Increaaes 
aaaociated with raglona of phaae tranaforaationa (AT(P.T,)) and thermal 
boundary layera (AT(B.L.))> reapectlvely [cf; Vetfhoogen, 1965; McKenzie, 
et al. 1974]. We conatraln first the left hand side of [1] and then 
evaluate each term on the right hand aide in succession. 

We presume throughout this discussion that the lower man, tie is 
convectlng, as is strongly indicated by the fact that it is very nearly 
adiabatic throughout [e.g. Jordan and Anderson, 1974; Dzlewonski, et al . 

1975; see also Stacey, 1977a]. Since an adiabatic gradient la not a 
steady solution to the conduction equation, it can exist only under very 
special circumstances and the most likely Interpretation i.s that the 
planetary interior is adiabatic due to sufficiently vigorous convection. 

The temperature at the top of the outer core is poorly constrained, however 
we estimate that a minimum plausible value is about 2800 K. We have estimated 
both liquid! and solldl for the Fe-S system under core conditions based 
on extrapolations of Usselman's [1975 a,b] high-pressure data for the 
eutectic and Liu's [ 1975] ‘ analysis of the melting curve of iron, corrected 
upward slightly to conform with the nonzero shear modulus observed in iron 
at 180 GPa under Shock [Al'tshuler, et al . 1971]. Taking Ahrens' [1979] 
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bound* on Che maximum sulfur concent allowed in Che core and assuming 
that Che core is at a temperature at least as high as the mean between 
eutectic and liquidus» we arrive at a lowest estimate of about 2800K 
(under these conditions Che outer core would already be a slurry con** 
slating of at least S-IOZ crystals according to our extrapolations). 

Our bound of 2800K is substantially less Chan most estimates of the temperacure 
in the core [cf. Jacobs, 1975], and is less than or comparable to previously 
suggested lower bounds [e.g*, Birch, 1972; Usselman, 1975b]. Furthermore, 
we have only considered sulfur as the light element in Che outer core: 
sulfur is likely to decrease the melting temperature of iron significantly 
more than other plausible candidates [e.g., Brett, 1976], in particular 
oxygen which is currently considered a strong candidate [Ringwood, 1978]. 

In this case a minimum estimate for T(2900) . would be above 3000-3200K, 
taking the zero pressure phase equilibria [Darken and Gurry, 1953] and 
assuming that the effect of pressure in the Fe-0 system can be scaled to 
that in the Fe-^S system [see also Ringwood, 1978]. Thus a more plausible 
value for T(2900) might lie between about 3200 and 3500K, in accord with 
recent analyses [e.g., Stacey, 1972, 1977b; Verhoogen, 1973; Usselman, 1975]. 

The temperature at the base of the transition zone can best be constrained 
petrologlcally, yielding estimates of about 1900-2000K for T(700) [Ahrens, 

1973; Aklffloto, et al . 1976]. For example, taking the transition from olivine 
to 8>splnel as beginning at about 390 km depth [Burdick and Helmberger, 

1978] provides a temperature "fixed-point” of about 1700K [e.g., Akimoto, 
et al . 1976]. This implies a temperature of about 2000K at the base of the 
transition zone, for plausible values of thermal expansion, specific heat 
and Clapeyron slope (cf. Verhoogen, 1965; Mao» et al . 1969; Schubert, et al . 
1975; Navrotsky, et al . 1979], Because of the sensitivity of the 


estimated temperatures on the detailed phase equilibria as well as 
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the uncerUlntlM in our eatlmatee oi' the thenw)- 
cheaical pro^iertleii In the transition xone, our values for T(700) are 
aoMwhat uncertain* perhaps by as much as 100-20GK. However a recent* 

In situ determination cf the ollvlne-splnel phase boundary [Aklmoto* et al. 

1977] suggests that our values of temperature may* If anything* be high. 
Furthermore* Independent evaluations of the temperature at the top of 
the lower mantle are very compatible with our result [e.g.* Banks* 1969* 

Graham and Dobrzykowskl* 1976; Watt and O'Connell, 1978]. Thus* we 
estimate the difference In temperature between the top of the lower 
mantle and the outer core to be no smaller than BOOK* and more likely of 
the order of 1300K. 

If we assume* initially* that the lower mantle is homogeneous (In 

phase and composition) the bracketed term In [1] should yield the rise 

In temperature from top to bottom of the lower mantle. Following Jeffreys 

3 

[1930]* dynamic stability requires a Rayleigh number less than about 10 * 

and this provides an estimate of the superadlabatlc temperature Increase through 

the lower mantle of IK. Here* the constraints currently available 

for the viscosity [Peltier and Andrews* 1976; see Table 1 for representative 

parameter estlmates]rhave been used and this estimate is only accurate to 

about one order of magnitude. In any case* the superadlabatlc gradient 

Is negligibly small. 

The adiabatic temperature gradient is given by the GrUnelsen parameter* 

* as a function of density (p). We have determined y for the 
S 

lower mantle from experimental data on MgO and ^^^2 throughout the pressure 
range of Interest. Shock-wave data on porous (or anorphous) and nonporous 
ssnples yield values of y directly as a function of pressure or density 


_ / 3lnT 
" ^31np 
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[■•e Garter t ct al. 1971, for example]; for MgO and SIO 2 these values 

are quite similar as a function of pressure [data base: Wackerle, 1962; 

Al'tshuler, et al. 1965; Hart and Skidmore, 1965; Anderson, 1967; 

Jonee, a t al. 1968; Ahrens and Rosenberg, 1968; Carter, et al. 1971; 

Trunin, et al. 1970, 1971; Pddurets, et al. 1976). Figure 1 shows the 

values of y derived from these data for mixtures of the oxides corresponding 

to olivine and pyroxene stoichiometry. The differences between these two 

cases are not resolvable In thin study. Taking MgO and SIO 2 at high 

pressures to be representative of the earth's deep interior, this supplies 

the ratio of temperatures at the top and bottom of the lower mantle along 

an adlabat by applying selsroologlcally determined density profiles 

[e.g., Dzlewonskl, et al. 1975]. Consequently AT , the adiabatic 

s 

temperature Increase across the lower mantle, corresponds to about 460K 
"" AXnT 

(with an average y 0.88), and the bracketed term In [1] amounts 

to no more than about 500K. This Is approximately 300K less than even 
our limiting estimate of the left hand side of [1]. 

That our estimate of the Grllnelsen parameter Is similar to those 
resulting from Independent analyses is shown In Figure 1. The different 
estimates do not strongly affect ATg In that Increasing y from 0.88 to 1.0 
raises the ratio of temperatures at the top and bottom of the lower mantle 
by only 3%. However, our analysis has the advantage of being constrained 
by experimental data measured at the appropriate pressures. 

Superadlabatlc Gradients 

In the discussion above we have estimated the temperature at a depth 
of 700 Km, a lower bound on the temperature at the core-mantle boundary 
and the GrUnelsen parameter, Which specifies the adiabatic gradient between 
these two points. Our results are Illustrated In Figure 2 and It Is apparent 
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that a alngla adiabat doaa not join the two tenperature eatiaatas. There- 
fore one or more reglona of auperadlabatlc gradients are likely to exist 
between 700 ka depth and the Core-mantle boundary. The required super- 
adlabatlc temperature Increase Is so large as to further validate our 
Initial assumption that convection Is Inevitable. We now explore how 
such superadlabatlc gradients arise In a convectlng mantle. 

If regions of phase transformation exist within the convectlng 
lower mantle » they could contribute superadlabatlc temperature Increases^ 
depending on the enthalpies of reaction Involved. The occurence of phase 
transformations has Indeed been Inferred from several selsmologlcal studies 
[a*g. Kanamorly 1967; Johnson, 1969; Vinnik and Nlkolayev, 1970; Wright 
and Cleary, 1972; Butler and Anderson, 1978] and, although the evidence Is 
still open to question [e.g., Wiggins, et al. 1973]^, as many as five seismic 
discontinuities have been Identified [Johnson, 1969]. Whether these 
transitions are univarlant or polyvariant, the additional change In 
temperature associated %d.th each transformation Is very nearly [Varhoogen, 
1965] 


AT(P.T.) (~) AV, 

Where T Is the temperature, C an appropriate (averaged) specific heat, 

^ the Clapeyron slope and AV the volume change associated with the 
transformation. The Inferred lower mantle discontinuities are known to 
be small (velocity jumps less than about IZ) and we take AV '^'2Z as a 
generous upper limit. For typical values of the Clapeyron slope (2-3 MPa 
we find AT(P.T.) no larger than 20-30K at each transition, and hence 
the total contribution from phase transformations to the temperature change 
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throughout: the lower aentle amouata to no more than about 100~130K, aa 
ahown In Figure 2. Of courae, if amaller (or negative) Clapeyron slopes 
are considered for any of the Inferred transitions (as might be appropriate 
for transformations involving disproportionation) the total AT(P.T.) is 
aiSaller than our calculated values (or even negative) . Clearly* the 
effect of phase transformations la insufficient to reconcile our two 
temperature estimates and we are therefore forced to consider further 
sources for auperadlabatlc temperature increases within the lower mantle. 

An obvious place for non-adiabatic gradients would be immediately 
above the core-mantle boundary if a significant fraction of the earth's 
total heat flow is supplied to the mantle from the core. Heat flowing 
out of the core will result in a thermal boundary layer at the base of 
the mantle and we need an appropriate convection model to estimate the 
characteristic thickness and temperature change across this boundary 
layer. One complication la that the required model will have to take 
Into account the effects of compressibility and spherical geometry when 
applied to such a deep layer but to date these two properties have only 
been studied separately* and in the case of spherical geometry over a 
very limited parameter range. Nevertheless the existing results are 
sufficient to evaluate the way in which the boundary layers from in- 
compressible* plane layer convection models will be modified. 

We will first consider the boundary layers that are found in plane 

layer calculations which assume negligible compressibility* uniform 

-2 -2 

viscosity and a surface heat flux of 5.85 X 10 Wm , which is close to 
ths average measured terrestrial heat flow. The material properties of 
the convecting material are given In Table 1. The resulting horizontally 
averaged temperature for three models of heating are given in Figure 3. 
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Wh«n all of the heat la supplied from belov* the averaged temperature 
structure coneista of two equal boundary layers joined by an isothenaal 
interior region. The characteristic boundary layer thickness 6 is about 
100 km and the aesociated temperature change AT(B.L. ) is approximately 
300K. 

A simple boundary layer analysis [see McKenzie, et al. 1974] can be 
used to determine the dependence of the boundary layer properties on the 
heat flux and material properties of the layer: 



where the symbols and their typical value for mantle materials are given 
in Table 1. Eqiuitlons [2] and [3] are consistent with the idea that the 
thermal boundary transmits the flux F conductlvely: 






F 


An Important point regarding the boundary layer properties is that 
uncertainties in the flux across a given boundary layer affect the estimated 
temperature change much more than the boundary layer thickness. This 
effect can be seen in Figure 3, which Includes a case in which half the 
heating is from below and half from within (by internal heat sources). 

The thickness of the lower boundary layer is not much different than in the 
case fully heated from below, but the temperature change is substantially 
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reduced. The fully Internelly heeted case (l.e., no flux from che core) 

Is also shown In Figure 3 to Illustrate the absence of any lower boundary 
layer under such conditions. 

The next question is how does spherical geometry affect the results 
given above? A useful calculation relevant to this question is by 
Young [1974] who considered finite amplitude convection between isothermal 
shells. Figure 4 compares the horizontally averaged thermal structure of 
convection In a spherical shell to that of convection between plane 
Isothermal boundaries. The spherical case has a relatively low Rayleigh 
number of about five times critical and a Prandtl number of five, which 
means that Inertial terms In the momentum equation are not negligibly 
small. The plane layer case has a somewhat larger Rayleigh number of about 
15 times critical and the Prandtl nisaber is assumed to be Infinitely large. 
The two calculations thus differ In more than geometry alone but the major 
effect on the thermal boundary layer structure does appear to be simply 
related to the geometric differences. The plane layer case has symmetric 
boundary layers. In the spherical case both the Inner and the outer 
boundaries transmit the same total amount of heat, but since they differ 
In surface area by a factor of about three, the heat flux per unit area is 
three times greater across the lower boundary. As one would expect from 
relations [2] and [3] this difference In heat flux has little effect on 
the boundary layer thickness but results In an almost three times greater 
temperature drop across the lower boundary layer. As a first approximation, 
the effect of spherical geometry can be taken into account by simply using 
the local heat flux for estimating boundary layer properties. 

The effect of compressibility is illustrated in Figure 5. The in- 
compressible temperature profile shown for reference Is the same as In 
Figure 3 but now using a specific depth extent of 700 Km. The thermal 
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structure la taken froB a recant atudy of convection In a coapreaalble 
■aterlal by Jarvla and McKenzie [1979]. The Inportance of compreanlblllty 
In auch a calculation la aeaaured by the ratio of the layer depth d to 
the thenaal acale height (H^ “ C^/ga). Thla acale height la the 
vertical dlatance a parcel of material muat be dlaplaced for Ita temperature 
to change by a factor e due to adiabatic ccnpreaalon. The caae used to 
Illustrate the effect of compressibility on the average thermal structure 
aasuses d/H^ - 0.5| which results In a reasonable large effect. Comparing 
Che two calculations one can see chat the Interior temperature la no 
longer constant In the compressible case, and In fact, the Interior 
temperature gradient Is found to be almost exactly adiabatic as we have 
taken It to be In our earlier discussion. The boundary layers are not 
vety much affected by the compressibility because they are small compared 
to When considering the mantle, we argue that compressibility will 
not greatly affect the thermal boundary layers because 6/H^ «l(6'v 100 Km, 

5000 Km) but the interior temperature gradient Joining the boundary 
layers must be considered to be adiabatic. 

Guided by the discussion on the effects of sphericity, we will 
assume that 1/4 of the total terrestlal heat flow arises originally In 
Che core. The heat flow per unit area out of the core is then almost 
exactly equal to the surface heat flow and we expect, based on the 
arguments given above, that the boundary layer at the base of the mantle 
will be about 100 km thick and with a teiDperature change of approximately 
300K. A smaller flux from the core would result In a smaller boundary 
layer, both in thickness and In superadlabatlc temperature gain 
(equations [2] and [3]). Similarly, because we have ignored the temperature- 
dependence of viscosity within the thermal boundary layer, our values of 
6 and AT may be somewhat large. However, Caking our best estimate of the 
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trapcracure at the top of the mantle, combined with our adiabatic gradient 
(± contributiona from phase transformations), we can satisfy estimates 
of the minimum core temperature once we include a thermal boundary layer at 
the bottom of the mantle. 

Independently, there is considerable evidence from seismology [see 
Cleary, 1974, for example] that the lowermost mantle [the D" region, 
extending to about 100-200 km above the core-mantle boundary: Bullen, 

1949] is distinct from the overlying mantle and does not behave as an 
adiabatic and homogeneous region. The seismologlcal anomalies corresponding 
to this region may therefore reflect the properties of a thermal boundary 
layer. We note that our estimated average temperature gradient through this 
layer ('v 3K km"^) is probably smaller by a factor of 2 to 4 than the critical 
gradient required to make the gradients in density or seismic parameter 
vanish; however dispersion may significantly affect the selsmologlcally 
observed velocities. Allowing the values quoted in Table 1 to Vary within 
reasonable limits does not significantly change our estimate of the thermal 
boundary layer thickness, therefore if the seismologlcal data require an 
anomalous layer thicker than about 100 km at the base of the mantle, this 
should be taken as evidence that something more than a single thermal 
boundary layer is present. 

Chemical Composition 

We have presented our best estimates of Che temperature field through 
the lower mantle and have found Chat a solution can be derived which is 
in agreement wlCh our minimum bounds on the core CemperaCure. We are 
required, however, to Inserc a slgnificanc chermal boundary layer at the 
core-mantle incerface (Figure 6a). This conclusion is not affected by 
the presence (or lack thereof) of phase transformations in the lower 
manCle because we found ChaC such CranslCions have a minor if noC negligible 
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•ffact on the pverell geothcrmil gradient. Ue note that the thenal 
boundary layer near the earth* a aurtace Involvea a nuch larger At than 
that at the baae of the oantle due largely to the teaperature dependence 
of the effective viscosity. 

Although this aodel is acceptable* it barely satisfies our bounds and 
we find it a rather extreme solution. If the lower mantle contains only 
one boundary layer* we estimate a temperature at the top of the core which 
is very low: about 2800 - 3000K. Such temperatures can not be considered 
unreasonable* however they are significantly lower than recent estimates 
[e.g.* Verhoogen* 1973; Stacey* 1977b; Bukowinski* 1977] and especially 
previous results ranging between 4000 and 5000K [cf. Jacobs, 1975]. 

One possible factor might come from our having underestimated the heat 

3/4 

flux from the core. Since AT(B.L,) is proportional to F (equation [3])* 
a larger flux might allow a large enough temperature Jump across the boundary 
layer to yield more plausible core temperatures. If the total flux from 
the core is comparable to the total heat flux through the earth's surface 
then our mantle temperatures are essentially compatible with our best 
estimates for the core. Such a large flow of heat from the core has* however, 
been generally considered unreasonable [e.g.* Stacey* 1972; Verhoogen* 1973]: 
our value of 1/4 of the total output at the surface represents the maximum for 
steady state. This Is within a factor of two of the estimated minimum [Stacey* 
1972* 1977a] and is probably a fairly good estimate. 

Thus, we find that if the lower mantle contains only one boundary layer 
at its base* the core must be at a significantly lower temperature than has 
previously been considered (essentially consistent with Tolland's [1974] 
conclusions). Alternatively* if a temperature of about 2900-*3000K Is 
considered to be too low for the outer core* the lower mantle must contain 
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oth«r boundary layers In addition to that at the core-nuintle boundary. 

He conclude, therefore, that a second boundary layer within Che lower 
■ancle Is suggested by our results up to this point. 

From our analysis we can not constrain at which depth such a 
boundary layer might be found, although we can rule out the possibility Chat 
there could be many such boundary layers within the lower mantle since the 
resultant temperatures In the core would most likely be well above Chi 
melting point of Iron at the Inner core ** outer core boundary. 

We consider two depths, as shown In Figures 6b and c, as the most 
plausible locations for such a boundary layer based on selsmologlcal evi- 
dence: either at the transition zone or near the base of the mantle. The 
lowermost tsantle, for eiumple. Is known to be anomalous and It has been 
Inferred to contain significant Inhpmogeneltles [e.g., Phlnney and 
Alexander, 1969; Davies and Sheppard, 1972; Kanasewlch, et al. 1973; 

Haddon and Cleary, 1974; Doornbos, 1976; van den berg, et al. 1978]. 

Thus, If layer D" Is chemically distinct from the overlying mantle It 
could contain a boundary layer both at Its cop and bottom, and the 
temperature could Increase by about 900K within the lowermost 200-300 km 
of the mantle, easily satisfying currently accepted temperatures for the 
core (Figure 6b). Such a model has. In principal, already been proposed 
[e.g., Sacks and Beach, 1974; Sn^k'^ and Sacks, 1976] and It has been 
shown to be compatible with some of the selsmologlcal data. However, 
numerous difficulties remain [e.g., Husebye, et al. 1976; Miller, et al. 
1977] and this model must be considered speculative at present. An attempt 
to explain the observations from the lowermost mantle, such as the velocity 
anomalies and apparent scattering. In terms of a slnq>le (single) boundary 
layer could provide an Important selsmologlcal test on Che constitution 


of Che lower mantle. 
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Th« rwuloing pottlblllty !■ to place a charaal boundary layer at Che 
cop o]E Che lower manCle (aaaoclaced wlCh Che CranalClon xone) aa veil aa 
ac, lea baee (Figure 6c) . Again, Chla poaalblllCy can noc be crlclcally 
aaaeaaed for lack of dacri, however ic haa laporcanc Inpllcaclona In chac 
IC precludea whole nanCle convecclon. Such a boundary layer would aoaC 
likely be aaaoclaced wlch a change In chenlcal conpoalClon froai upper co 
lower nancle [cf. RlchCer and Johnaon, 1974] and la cherefore conalacenc 
wlCh auggeaClona Co Chla cffecc. Indeed, If Che IndependenC arguaenca In 
aupporc of a change In coapoalclon acroaa Che CranalClon cone are correcC 
[e.g.» Anderaon, 1968, 1977; Preaa, 1968; Anderaon and Jordan, 1970; 

Anderaon, eC al. 1972; Burdick and Anderson, 1975; Liu, 1979], Chen a 
boundary layer would be expecCed ac Che Cop of Che CranalClon zone and 
Ceatperacurea Chrough Che lower nanCle would be aome 500-700 K higher Chan 
we have ahown In Figure 2. WlCh regarda Co experimenCal, high preaaure pltaae 
equilibria, Che aicuaclon aC preaenC la thac no adequate explanaClon haa 
apparenCly been found for Che large veloclcy and denalcy jump ac Che base 
of Che CranalClon zone (e.g.. Hart, e t al. 1977; Burdick and Helmberger, 

1978] In terma of polymorphic phase tranaformaclona alone [Anderson, 1976; 
Liu, 1979], and therefore a change In bulk chemical composition Is strongly 
Implicated for the transition zone. Similarly, recent Isotopic and trace 
element data can be Interpreted to suggest the existence of a separate, 
relatively undepleted geochemical reservoir underlying a more depleted 
upper mantle reservoir [DePaolo and Waaserburg, 1976; O'Nlona, et al. 

1978* 1979]. 

Conclusions 

The central point of Chla study la Chat there Is good evidence for 
regions of non-adlabaclc gradients within Che lower mantle and these ace 
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•Oft likfly ton£lnid to thtmal boundary laydrs. W« hfvf ovaluatfd 
thf adlfbatio tharwil gradlfnt through th« lower mintlo by ualng high 
priffur* •xporinentfl data on tha GrUnalsen paraaatar and have found that 
thla gradient if cbapatlble with a lower bound on the temperature in the 
core only after allowing for a thermal boundary layer at the core-mantle 
boundary. Even with a thermal boundary layer at the core-mantle boundary 
we require a very low core temperature^ and therefore we dlfcueaed the 
poaaibility of further non-adiabatic regiona within the lower mantle. 

Our analyais teliea upon having properly estimated the temperature at 
the top of the/mMtie and other constraints on the nature of the transition 
xone from high-pressure experimental petrology would be very useful in this 
regard. However our geo therm can be arrived at by oeveral lines of reasoning 
and it is therefore considered to be adequately constrained. 

We can not resolve the effects of possible phase transformations 
within the lower mantle since they appear to have only minor consequences 
for its thermal state, within the scope of our study. The thermal boundary 
layer properties depend on the thermal flux assumed from the core, and our 
conclusions regarding the need for more than one thermal boundary layer within 
the lower mantle, may not be valid if a significantly larger flux is allowed 
from the core than has previously been considered acceptable. Alternatively, 
we find that either the temperature in the core must be lower than has been 
considered plausible, or there must be a barrier to convection within the 
lower mantle which results in additional thermal boundary layers. We find 
this latter possibility more plausible, however we can not constrain the 
depth at which such a barrier may be found. Guided by selsmologlcal evidence, 
we suggest the existence of a transition in chemical composition either at 
the top of the lower mantle or toward its base (D" region) as being the 
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aoft likely cause o£ extra thermal boundary layers. ThlSf :huni provider 
ua vlch a large enough Increase In teaperature across the lower taantle to 
allow plausible core teaperaturea (about 3800K) . If a coaplex zone with 
aultlple theraal boundary layers and .separate convection within it la 
required to explain the observed (anaaaloua) properties of the D" region 
then we eatlaate an Increase in teaperature of about lOOOK within f.he 
loweraost 200-300 ka of the aantle. Alternatively, the transition zone 
is a likely candidate for a region across which the bulk chealcal coaposition 
changes; this could readily explain a theraal boundary layer located there. 
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Table 1 

Lover Mantle Properties 


Density: 

P 


4 

to 

S.5 Mg 

m 

-3 

(1) 

Acceleration of Gravity: 

g 


10 ms“^ 


• 

(1) 

Specific Heat at Constant Pressure! 


'U 

1.2 X 10^ 

J 

Kg“V^ 

<2) 

Kinematic Viscosity: 

V 


1 

to 

17 2 -1 

10 X 10 ' s 

<3) 

Thermal Dlffuslvlty: 

K 

% 

1 

to 

2 X lO" 

•6 

2 -1 
m s 

(4) 

Coefficient of Thermal Expansion: 

a 


1 

to 

2 X lO’ 

•5 

K-' 

(5) 

Neat Flux From Core: 

F 

m 

5. 

85 

X 10*"^ 

W 


(6) 


Sources 

(1) Measured 

(2) Dulong-Petlt value 

(3) Peltier and Andrews [1976] 

(4) Schatz and Simoons [1972]; Kleffer [1976] 

(5) This study; a - ^ ; also, Stacey [1977b] 

(6) Assumed equivalent to the flux through the earth's surface: 


see text 
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Figure Captiomi 

Figure 1 GrUneieen parameter through the lower mantle based on 
experimental high pressure data for MgO and SIO 2 (this 
study) compared with the theoretical formulations of 
Stacey [1977b], and Brennan and Stacey [1979]. 

0 and P refer to mixtures of the oxides corresponding to 
Olivine and pyroxene stoichiometry » respectively. 

Figure 2 The temperature profile for a homogeneous and adiabatic 

lower mantle (solid line) is compared with profiles including 
the effects of possible phase transformations (dashed line) 
and a thermal boundary layer at the core-mantle boundary 
(dotted). Fixed reference points in the transition zone 
are shown* along with our best estimate and lowest plausible 
values of temperature in the core. 

Figure 3 Horizontally averaged temperature as a function of depth 

from nianerlcal simulations of convection In an incompressible 

fluid layer. [McKenzie* Roberts and Weiss* 197A]. The flux 

-2 -2 

at the upper surface Is 5.85 X 10 WIm and relevant material 
properties are given In Table 1. 

a) All heating from below 

b) All heating internal 

c) Half heating from below* half internally heated 

The boundary layer thickness S is measured from the boundary 
to the first point where dT/dz ■■ 0. 
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riBuf 4 Horliontally avaragad taaparatura froa nuMrlcal cmlculationa 
of convactlon batwaan laotharmal boundariaa. 

a) PXana layar': Ra - 10^, Pr - 

b) Spharical ahall: Ra ■ 4 X 10^, Pr ■ 5 and ratio of 
Innar radlua to outar radiua of 0.6 (axially ayaaatric 
■oda, takan fron Young, 1974] 

3 

Ra ■ , Pr ■ v/k and AT ■ T 2 - Tj^ 

The vertical axla Is depth in Cractlona of d, tha total 
layer depth. 

Figure 5 Horizontally averaged temperature for both Incompraaalbla 
and compreaaible convection. 

a) Sane aa in Figure 3a but with depth of 700 Km 

b) Conmresaibla caaa [Jarvia and McKenzie, 1978: 
case B 10] 

Assumes depth " 700 Kn, d/H^ ■ 0.5 and flux of 5 X 10**^ Me”^ 
from below. 

Figure 6 Temperature profiles through the earth's mantle: three 

acceptable models. Thermal boundary layers (B.L.) are indicated 
in the upper and lower mantle, while dynamic barriers associated 
with possible chemical transitions in the lower mantle are shown 
as dashed lines. The preferred core temperature and its plausible 
lower bound are also shown. Note the fundamental asymmetry of 
the profiles caused by the large decrease in temperature through 
the lithosphere: an effect which has no counterpart at depth 
within the earth and which is due to temperature-dependent 
rheology. 
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riKuf 6 (continued) 

(a) Adlabac through lower mentle with one themal bounder/ 
layer at the core-nantle boundary and with the core at 
a relatively low teaperature. Taken froa Figure 2. 

(b) Adlabat with two theraal boundary layera near Che baae 
of the nantle correapondlng to a aelaaologically coaplax 
and cheaically diadnet D*’ region* 

(c) Adlabac with a theraal boundary layer at the top and the 
bottom of the lower aantle. In thla caae, the lower and 
upper aantle are aeparace chealcal and dynamic ayeteBa. 
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PLATE TECTONICS ON VENUS 
Don L* Anderson 

Selsmologlcsl Laboratory, California Institute of Technology, 
Pasadena, California 91125 


Abstract 

The high surface temperature of Venus implies a permanently buoyant 
lithosphere and a thick basaltic crust. Terrestrial style tectonics 
with deep subductlon and crustal recycling is not possible. 
Overthickened basaltic crust partially melts instead of converting to 
ecloglte. Because mantle magmas do not have convenient acesss to the 
surface the ^^Ar abundance in the atmosphere should be low. Venus 
provides an analog to Archean tectonics on the Earth. 

Introduction 

The surface temperature of Venus is about 450 K warmer than the 
surface of the Earth. This affects the buoyancy, thermal expansion, 
thermal conductivity and, hence, the thermal evolution and ultimate fate 
of the lithosphere. The pressure In the upper mantle of Venus Is at 
least 12% less than at equivalent depths In the Earth's upper mantle. 
Thus, the depth of melting, the locations of upper mantle phase changes 
and the viscosity of the upper mantle will be different for the two 
planets. The buoyancy and thermal properties of the lithosphere control 
the style of plate tectonics and the associated time and length scales. 



Mott dltcuttlont! of the comparative tectonics between Earth and Venus 
sddrett only the differences in viscosity. 

The oceanic lithosphere of the Earth cools as it ages and it 
eventually becomes denser than the underlying mantle. This instability 
develops after about 40 m.yt (Oxburgh and Paramentler, 1977). Once 

cooling reaches a depth of 50 km the garnet pyroxenite or eclogite 
transformations may make a substantial contribution to the negative 
buoyancy of the lithosphere. Continents and oceanic plateaus resist 
subductlon because of their thick low-density crusts. The fate of the 
terrestrial lithosphere, therefore, depends on both chemistry and 
temperature. The purpose of this paper is to Investigate the 

Implications for Venus tectonics of its high surface temperature. We 
will conclude that the surface thermal boundary layer on Venus Is 
permanently buoyant and that the reversible part of mantle convection 
occurs below about 100 km. 

The terms lithosphere and thermal boundary layer will be used for 

the cold outer layer of a planet. Neither term is strictly correct 

since the former has a strength connotation and the density of the 
latter is usually assummed to depend on temperature alone. The 
thickness of the mechanical or rheological lithosphere depends on 
temperature and, possibly, composition and is proportional to the 
thickness of the thermal lithosphere that is discussed here. 

Lithospheric cooling 

The average temperature of the oceanic lithosphere decreases about 
660*C as it ages. The thickness of the conductively cooled thermal 
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boundary layer Increases at a rate controlled by the thermal 

conductivity and the difference between the surface and Interior 
temperatures. Since the crustal and harzburgite portions of the 

boundary layer are both less dense than the underlying mantle a 

gravitational instability can only occur after the conductive cooling 
penetrates for a sufficient distance, ~3Q km, into the denser upper 
mantle portion of the boundary layer. In some petrological models the 
lower lithosphere is pyroxenltlc. Upon cooling, plagloclase and spinel 
pyroxenite convert to garnet pyroxenite or ecloglte with a substantial 
increase in density. These phase changes also require cooling to depths 
in excess of 30 kilometers. 

The equilibrium thickness, 6, of a conductive boundary layer is 

6 • KAT/Q 

where K is the thermal conductivity, AT/ 5 is the average thermal 
gradient and Q is the mantle heat flow. For K of 7 x 10 cal/cm sec*C, 
AT of 1300®C and Q of 0.6 ycal/cn^sec, values appropriate for the 

terrestrial oceanic lithosphere, 6 is about 150 km. For a lithosphere 
composed of 6 km of basalt, 24 km of hatzburglte and a lower lithosphere 
composed of undepleted mantle, the equilibrium boundary layer for AT ■ 
1300*’C is 1% denser than underlying mantle and therefore gravitationally 
unstable. Phase changes in the lower lithosphere may contribute further 
to the negative buoyancy. 

For a given mantle temperature and heat flow the equilibrium 
thickness of the boundary layer on Venus is reduced by two effects; K 
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Is about 30% lover at the higher teaperature (Schetz and SlazDona, 1972) 
and AT is reduced by 460*C^ This reduces 6 on Venus to 42 km and its 
density, assuroining similarity with oceanic lithosphere, is 2% lighter 
than the underlying mantle. This is a result of the smaller amount of 
cooling, the higher proportion of basalt and harzburgite in the boundary 
layer and the lack of significant cooling where it is most required, 
l.e. the lower lithosphere. Furthermore, phase transformations in the 
lower lithosphere cannot contribute at high temperature and low 
pressure. Basalt has about one-half the thermal conductivity of 
ultramafic rocks. This reduces the thickness of the boundary layer to 
about 27 km if it is mainly basaltic. 

The buoyancy of the Venus lithosphere, even for a relatively thin 
crust, is greater than the buoyancy of young oceanic lithosphere on 
Earth. Since a thicker crust for Venus is probable. It is certain that 
the surface thermal boundary layer for Venus is permanently buoyant and 
stable against subductlon. 

The temperature gradient in the conductive layer can be written 

dT/dz - (l/K)(Q - Az) <2) 

s. 

where Q is the surface heat flow and A is the heat production rate in 
s 

the layer. With basaltic conductivities and radioactivities (A ~ 10“^^ 
cal/cra"* sec) the temperature and thermal gradient at 20 km are 950"C and 
23"C/km respectively and the solidus of dry basalt will be exceeded at 
depths shallower than 40 km. Eclogite is not stable in the mantle of 
Venus at depths shallower than about 100 km (Anderson, I960). 
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Since crust cannot subduct and the garnet-rich asseablages are well 
below the boundary layer, all partial melt products of mantle 
dlffercntation that make their way Into the outer 100 km or so will 
remain there. The crust is therefore likely to be thi'^ker than either 
the oceanic crust or the average crustal thickness on Earth. It may be 
thicker than either the mechanical lithosphere or the thermal boundary 
layer. Ctust which is overthickened by compression, buckling, 
thrusting, or collision will melt at its base. The resulting plutonism 
will further increase the density contrast between crust and mantle and 
increase the radioactivity and incompatible trace element content of the 
parts of the crust so affected. 

The characteristic time for thickening of the boundary layer is 

T - (pCp/4 K)62 <3) 

which for p ■ 3.3 g/cm^, Cp « 0.25 cal/g*C and 6 * 150 km is ■^200 m.y. 
for the Earth. For Venus the characteristic time is 20 m.y. or less. 

Therefore, Venus has a thick, low-density and permanently buoyant 
crust-lithosphere that very quickly reaches thermal equilibrium. Upper 
mantle temperatures are high, resulting in low viscosities. 

The basalt-ecloglte transformation will occur below about 100 km. 
The driving mechanism for mantle convection may therefore be partly 
chemical. Basaltic melt rising in an upeurrent will freeze as eclogite 
if it remains below about 100 km, and provide the negative buoyancy 
required for overturn. The 670 km discontinuity which halts eclogite 
subduction on Earth (Anderson, 1979a, c) will occur at about 800 km in 
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V«nus giving a convactlng layer about 700 km thick* Thla ahould 

ganarata topographic and gravity anomallaa of about 1400 ka Vavalength* 

Equation (1) auggcate that plate tactonlca In the Archaan on Earth 

■ay have been similar to present day tectonics on Venus. The small i In 

early Earth history was primarily due to the high Q from radioactive 

decay, which was at least three times larger 4 x 10^ years ago. The 

thermal boundary layer thickens as Q decreases until eventually It 

becomes unstable and overturns, destroying the early geological record 

and setting the stage for the present style of tectonics (Anderson, 

1979a, b, 1980). This overturn event would allow much of the accumulated 
AO 

Ar In the upper mantle to escape and the subsequent steady-state 
ridge/trench style of tectonics allows continuous outgasslng of the 
terrestrial mantle. 

Early Evolution of Earth and Venus 

The Initial evolution of Earth and Venus were likely quite similar. 
The high accretlonal energies would result In partial melting of the 
mantle and upward transport of a plcrltlc melt. Crystallization of this 
upper mantle melt layer, or magma ocean, yields a thin plagloclase-rich 
crust and a deep ecloglte cumulate layer. An ollvlne-orthopyroxene 
ciimulate layer may form at Intermediate depths. This kind of scenerlo 
Is well documented for the Moon except that pressures are too low for 
extensive ecloglte fractionation. Ecloglte fractionation removes AI 2 O 3 
from the melt and reduces the thickness of an early anorthositlc crust. 
Ecloglte cumulates sink no deeper than 670 km on Earth because of the 
Intervention of the llmenlte and perovsklte phase changes In a 
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pcrldotltic iuntl« (Andefson^ 1979c). An ecloglCe cumulate lA’/tc wil)< 
ba daapar and amallar In Vanua because of the effect of temperature and 
praaaure on the phese boundaries and the limited stability field of 
aeloglte In the upper mantle. Nevertheless, the early tactonlca, 
gaocheroical dlfferentatlon and outgasslng of the two planets wars likely 
to have been similar. 

If the outer layer becomes unstable as the planet cools, the 
situation changes. Subducted slabs affect both the flow and thermal 
regime of the upper mantle. The sinking of the cold boundary layer deep 
Into the upper mantle drags cold Isotherms to depth and allows hot 
replacement material to rise at ridges. This leads to relatively large 
lateral and vertical temperature differences In the upper mantle and 
rapid convection in which the surface boundary layer participates. 

The mantle of the Earth contains approximately 15 ~ 20 % of a basaltic 
component. If the mantle were well differentiated this would give a 
crustal layer 350**450 km In thickness. Only the top 50 km would be 
buoyant, because of the ecloglte transformation, but this Is still more 
than 4 times the present average crustal thickness. The lunar and 
Martian crusts are much closer to being the appropriate thickness for a 
well differentiated planet in spite of the fact that accretlonal heating 
must have been less for these bodies. The difference, of course, is due 
to Che continuous removal of basaltic crust from the surface Implying a 
long-term storage reservoir at depth. The size of this reservoir has 
been estimated to represent about 20% of the mantle (Anderson, 1980). 




Samples thought to be representative of the upper mantle are 
depleted in AI 2 O 3 , CaO and SIO 2 compared to cosmochemlcal estimates of 
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■anfclt conposlClon. Thty art alto dapletad in Craca alamanta which are 
retained by the eclogite eineralsi garnet and clonopyroxene. The degree 
of depletion implies an eclogite layer about 450 km thick which is 
probably deeper than 200 kot 

The ^^Ar abundance for Venus is about an order of magnitude less 
than for the Earth (Pollack and Black, 1979). This suggests that late 
outgassing has been less efficient for Venus than for the Earth in spite 
of the higher surface and upper mantle temperatures. This is easily 
understood with the present model. On Venus, mantle melts and their 
volatiles mainly have access to the surface through relatively transient 
rifts. Because of the conservation of near surface material these rifts 
can only form when the crust compresses or thickens somewhere else on 
the planet. The amount of early outgassing on the two planets may have 
been similiar but Venus did not experience the early lithospheric 
overturn event or the outgassing asssoclated with continuous crustal 
renewal at ridge axes. 

Conclusions 

The high surface temperature of Venus has several important 
tectonic Implications. The most significant Is the small amount of 
cooling that the lithosphere experiences before it reaches thermal 
equilibrium. This increases its buoyancy and long-term stability at the 
surface when compared with the terrestrial oceanic lithosphere. Since 
crust cannot be recycled the mantle the thickness Of the basaltic 
crust on Venus is much greater than on Earth. The combination of a 
thick crust and high temperatures decreases the thermal conductivity of 
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the lithosphere. This means that for a given mantle heat flow the 
temperature gradient In the lithosphere Is greater than would be the 
case for a cold lithosphere with a thin crust. This leads to high upper 
mantle temperatures, low viscosities and the possibility of partial 
melting at relatively shallow depths. 
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Abstract 


Trace element c.oni;entratlQn patterns of continental and ocean 
island basalts and of mid**ocean ridge basalts are complataentary . The 
ralativc sizes of the source regions for these fundamentally different 
basalt types can be estimated from the trace clement 
enrichment/depletion patterns. Their combined volume occupies most of 
the mantle above Che 670 km discontinuity. The separate source regions 
are the result of early mantle differentiation and crystal fractionation 
from Che resulting melt. The MORB source evolved from an eclogite 
cumulate which lost its late stage enriched fluids at various times to 
Che shallower mantle and continental crust. The MORB source is 
primarily garnet and clinopyroxene , while the continental and ocean 
Island basalt source is a garnet peridotite chat has experienced 
secondary enrichment. These relationships are consistent with Che 
evolution of a terrestrial magma ocean. 



And«rion 3 


Introduction 


Hotspots* or plumss, hsve not yet been fitted satisfactorily Into 
either the tectonic or geochemical framework for the evolution of the 
mantle. From the point of view of mantle processes and the chemical 
evolution of the mantle the alkali volcanlsm associated with hotspots Is 
a very Important process even If the abundance of such rocks Is small 
compared to the volume of abyssal tholelltes (1). The processes of 
creation and subductlon of oceanic crust and lithosphere clearly account 
for most of the mass transport into and out of the upper mantle and the 
study of mid-ocean ridge basalts (MORB) has placed Important constraints 
on the part of the mantle that Is providing these magmas. Trace element 
and Isotopic considerations show* however* that the source region for 
MORB cannot be representative of the average composition of the mantle 
nor can It have existed as a separate entity with its present 
characteristics for the full age of the Earth. The composition, 
location and volume of this source region, which we designate MORBS, are 
still uncertain. It appears to be very large, homogeneous, global in 
extent, and to be the result of a previous differentiation or 
fractionation process that depleted It, relative to other source 
regions, of most of the incompatible elements. On the ocher hand, it is 
clearly not lacking in a basaltic component. The continental crust is 
complementary to MORB in both trace elements and isotopic ratios (2,3). 
Formation of the continents has presumably been at least partially 
responsible for the depletion of the so-called oceanic mantle. If the 
continental crust Is the only enriched reservoir* then mass balance 
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calculations Indicate that only about 30Z of the mantle has been 
processed, the rest remaining undifferentiated and primordial (3,4)* A 
large primitive reservoir would be difficult to reconcile with other 
evidence regarding the early thermal evolution of planetary Interiors. 
There is, however, evidence for an additional enriched reservoir in the 
mantle. Magmas from continental interiors, rift zones, oceanic islands, 
anomalous ridge segments and Island arcs, mantle xenollths from 
kimberlites and alkali basalts, and kimberlites themselves, all indicate 
the presence of a mantle reservoir with trace element concentrations 
that are also complementary to trace element concentrations in MORB. 
This reservoir apparently is also global and can provide magmas to a 
variety of tectonic environments including all types of plate boundaries 
as well as plate interiors. These magmas, of which alkali basalts are 
one example, are not nearly as voluminous as MORB but they indicate the 
presence of an enriched region or layer which may be substantial 
compared to the continental crust. We shall designate this as the 
hotspot or plume reservoir, or simply PLUME (S), and assume that it is 
the source of alkali basalts, nephelinites , melllltltes, basanltes, 
kimberlites and continental tholelites which it provides by varying 
degrees of partial melting (6,7). One of the diagnostic signatures of 
these basalts is the rare**Earth element (REE) pattern which Implies a 
source enriched by about 10 in LREE and 3 in the HREE, relative to 
chondrites (6,7,8). 

PLUME basalts also have high ratios of Rb/Sr, Rb/K, Ba/K, Nd/Sm, 
Ba/Nd, Ba/Sr, Nb/Zr, ®^Sr/®^Sr, and ^^^Nd/^^^Nd, and high concentrations 
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of H 2 O, CO^i Tl» K, Rb, Sr, Ba, La, Nd, Th, and U raXative to M0R6S. 
Thaaa charactariatics are shared by continental and ocean island basalts 
and are also evident, albeit diluted, In back-arc basins, Island arcs, 
and transitional or anomalous ridge tholelltes. Ocher characteristics 
that nay be common but which have been studied In only a few hotspot 
locations ar* high abundances of Cl, F, Br and the primordial rare gases 
(9). The process that led to the enrichment of PLIJME will have depleted 
other, presumably deeper, regions of the mantle. Tlie absolute and 
relative sizes of these enriched and depleted reservoirs bear directly 
on the problem of the evolution of Che mantle and, in particular, on Che 
question of whether there is a large primitive reservoir In Che mantle. 

Basalt source regions 

The continental crust is extremely enriched in the incompatible 
trace and minor elements. Continental tholelltes, basanltes, 
nephelenltes, alkali basalts, kimberlites and basalts from oceanic 
Islands are also enriched. The complement to the depleted mid-ocean 
ridge source region may, therefore, be much more voluminous than just 
Che continental crust. The trace element and Isotopic affinities of 
these rocks, regardless of their tectonic setting, suggests that they 
are all derived from a similar source region. This reservoir, PLUME, 
has previously been referred to as continental mantle, Che enriched 
source region or primitive mantle. Its location Is uncertain but it 
appears to be a global layer. Proposals for Its location Include the 
uppermost (6,10) and lowermost (11) mantle. In trace element ratios 
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•ueh «f Ba/Nd, Hd/Sa, B«/Sr» U/K and Rb/K there is a progressive 
increase fron ocean island and Island arc basalts, basalts from 
anomalous ridge segments, continental tholelites, alkali basalts to 
kimberlites* Mld-pcesn ridge basalts have much lower ratios and 
tholelites from back-arc basins are generally slightly higher chan MORB. 
Trace element concentrations vary systematically through the 
compositional spectrum olivine melilitite, olivine nephelinite, 
basanite, alkali basalt, and olivine tholeiite and can be explained by 
varying degrees of psrtlal melting ranging from 4 to 25% from a common 
source region which is enriched in the strongly incompatible elements 
(Ba, Sr, Th, U, LREE) by lOx chondtitic and in the moderately 
Incompatible elements (Tl, Zr, Hf, Y, HREE) by a factor of about 3 
(6-3). Mid-ocean ridge basalts are derived from a very different source 
which is depleted in Che large— ion lichophlle (LIL) elements* The 
complementary nature of the two source regions is possibly due to the 
migration from the MORB source region of a melt or fluid with 
incompatible element concentrations similar to those of kimberlite* 
This fluid depletes the MORB source region and enriches the 
complementary mantle reservoir* The upper mantle low-velocity zone 
(LVZ) is a likely repository of these volatile and trace element 
enriched fluids and, therefore, a possible source region for the 
enriched magmas which are characteristic of the plume or hotspot source* 
Kimberlites come from depths as great as 220 km, consistent with an 
origin near the bottom of the LVZ* The xenollths in kimberlites and 
alkali basalts represent mantle fragments from various shallower depths 
and also generally exhibit enrichment of the type inferred for the plume 
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•ouffc* rtglon. Th«r« li now a waalch of avidianea for upparaantla 
■acaiODUitlsin (12» 30) and anrlchaane avanta afc various elowa. Tha 
•ourea of tha anriching fluid has not baan isueh diacussssad. 

The coaplamantary nature of continental and oceanic tholalites is 
illustrated in tha canter part of Figure 1. Alkali! basalts, 
naphalinites, malllltitaa and basanitas have similar but mote enriched 
patterns than continental flood basalts, consistent with their 
derivation from the same source region by smaller dageas of partial 
melting. Continental and abyssal tholalites exhibit reflection symmetry 
about a line corresponding to a 6 to 7 fold anriehsMnt over average 
mantle concentrations. 

Continental and ocean island basalts are relatively depleted in 
those elements that are retained by garnet and clinopyroxene. The 
reverse is the case for abyssal tholelites. Since garnet and cpx are 
reduced or eliminated during the large amounts of partial melting 
usually inferred for the formation of tholelites (6), the Inverse ga-cpx 
pattern of plume basalts must be related to a prior history involving 
eclogite fractionation or invasion of PLUME by a fluid that was in 
equilibrium with a ga-cpx assemblage. 

The origin of island arc basalts is still controversial but they 
have PLUME affinities in such trace element ratios as Rb/K, Rb/Sr, Ba/K, 
Ba/Sr and Ba/Nd. These ratios are closer to ocean island basalts, 
alkali basalts and continental flood basalts than they are to raid-ocean 
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rldgt ba««\ta. Back»arc b4^.aln baaalCf are Incarmadiaea td MORfi and 
hotepot nagnaa, auggeatlng a mixture from a deeper MORB aource and an 
overlying enriched aource region. Ocean ialand tholeiitea alao appear 
to be mixturea both in trace element and iaotopic ratioa (2). The fact 
that both HORB and batilta with continental affinitiea are available at 
ridge enrl/OMMnta, island arca» oceanic islands and well developed 
continental rifts suggests a compositionally stratified mantle. The 
order of app«^arance of these basalt types is consistent with the 
shallower layer being enriched in the incompatible trace elements. In a 
stratified mantle the deeper layer should be denser and, therefore, more 
garnet rich. The. ratio of concentrations in MORB to those in 
continental tholeiites indicates that the MORB source region, although 
depleted in most of the incompatible trace elements, has selectively 
retained those that are most comfortable in the garnet and clinopyroxene 
lattices. TIUh suggests that MORBS is an eclogite or garnet pyroxenite 
cumulate. 

The composition and volumes of the two source regions 

The concentrations of some key incompatible trace elements in 
various magma types, normalized to average mantle concentrations, are 
shown in Figure 1. Relative to the average mantle, continental 
tholeiites, alkali basalts, kimberlites and the continental crust are 
all most enriched in K, Rb, Ba, La and U. These are the elements that 
are most discriminated against by the major mantle minerals, olivine, 
or'thopyroxene , clinopyroxene and garnet, i.e., they have the lowest 
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mln«ral/Mlt p«rtltldn coefficients* The otheir elenente ere eXso 
■trongly cejeeced by olivine and orthopyroxene, having partition 
coifficlenta of lees then 0.03 for these minerals (6). On the other 
hand, garnet and clinopyroxene have partition coefficients greater than 
about 0*1 for Y, Nd, Sm, Sr, and Yb* These are the least enriched 
elements in the magmas vith continental and ocean island affinities and 
Che most enriched in MOAB* It appears that garnet and clinopyroxene are 
abundant minerals in the MORE source region, and that the fluid which 
enriched PhUME was previously in equilibrium with a ga-cpx rich 
reservoir. 

The high concentration of Y and Yb in MORB relative to other 
elements and relative to concentrations in basalts from PLUME, Che 
enriched source region, is particularly significant. These elements 
have mineral/melc partition coefficients greater than unity for garnet. 
Mass balance calculations suggest tltat most of the terrestrial inventory 
of Y and Yb may be In MORBS* This can only be accomplished if garnet is 
a dominant phase and if most or all of Che mantle has experienced 
differentiation. Partial melting of primitive mantle concentrates Y and 
Yb Into the melt; crystallization of this melt at moderate presssure 
would concentrate these elements Into an ecloglte cumulate. 

The relative volumes of the two source regions can be inferred from 
Figure 1 if it is assumed that they, plus the continental crust, give an 
undifferentiated terrestrial pattern for the very Incompatible elements. 
Using 0.56Z for the mass of continental crust relative to the mantle (3) 
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w« compuce frhat th« maia of HORB ausc txc««d the thoLelltlc fraction of 
PLUME by a factor of 12 to 15. Continental tholelltea represent about 
20Z BMltlng of a perldotltlc source region (6) and MORB must be nixed 
with at least 15Z olivine to recover the conposltion of Its parent magaa 
(13,14). Tlie residual crystals In both cases are assuned to be depleted 
because of their very low partition coefficients (6,8). The MORB source 
region is, therefore, at least 2.B to 3.5 times the mass of PLUME and 
the enrichment of the combined source regions, relative to primitive 
mantle, is about 3.3 to 3.8. Note chat MORB Is depleted only in a 
relative sense. Compared to average mantle abundances it is enriched. 
This suggests that MORBS plus PLUME are complementary to the remainder 
of the mantle from which they have presumably been removed by partial 
melting. The enrichment factor implies whole mantle differentiation and 
a residual mantle about 2~l/2 times the size of the combined reservoirs. 

By assuming that the very incompatible elements have been entirely 
fractionated into Che two source regions it is possible to estimate 
their total siass. For example, if La is depleted in the lower mantle 
the combined source regions represent 26 to 30% of the mantle. Tills 
corresponds to a thickness of 560 to 640 km in the upper mantle, or a 
region extending upwards from the 670 km discontinuity to a depth of 27 
to 110 km. If we take the transition region, 220-670 km, to be the 
depleted source region (20), then PLUME is 130-160 km thick, about the 
average thlckneas of the low-velocity zone (LVZ). 


A similar calculation for the other strongly incompatible elements 
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<Kf Rbt U) gives s range of 21-38Z for ftORBS plus PLUME* These are 
upper bounds since It has been assumed that the whole mantle has been 
processed and depleted. This, however, is a good approKlmatlon if 
olivine and orthopyroxene are the main residual phases. It Is 
significant that the above estimates of the relative and absolute sizes 
of the principal mantle reservoirs correspond to the main subdivisions 
of tite upper mantle, the LVZ and the transition region* The 670 ka 
depth corresponds to a major seismic discontinuity and the maximum depth 
of earthquakes. 

Other elements can be used to estimate the sizes of the individual 
reservoirs* The partition coefficients of Yh and Y are such that they 
strongly prefer the garnet structure (15). Both are enriched in MORB 
relative to PLUME, suggesting a higher portion of garnet in the former. 
If the entire mantle complement of Y and Yb reside in MORBS, then this 
would represent lA-lIZ of the mantle or an upper mantle equivalent 
thickness of 300 to 364 km. Sc is partitioned into garnet and 
cllnopyroxene by a factor of 8 relative to olivine and ortliopyroxene. 
Using this ratio of enrichment for the MORB source relative to the rest 
Of the mantle it would constitute 26Z of the mantle. 

K, Rb, Ba, La and U have the highest relative concentrations in 
PLUME* These are the elements with the lowest partition coefficients 
for the major mantle ^nerals and those that are most likely to be 
concentrated into PLUME by melts or metasnmatic fluids* The PLUME 
abundances of these elements can account for the entire mantle Inventory 
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If PLUME is 9 to 19% of th« mantlo* 

Using the above estimates of the relative sizes of the two source 
regions and Ganapathy and Anders (37) values for terrestrial abundances, 
we can now estimate the total abundances of the trace elements In the 
various reservoirs# The continental crust contains more than 26* of the 
mantle plus crust Inventory of K, Rb, and Ba and less than 8% of Sr, Y, 
Zr, Nb, Sm, Yb, and Hf# The ''depleted" MORB source contains more than 
50% of the Earth's Inventory of Y, Zr, Nd, Sm, Yb, Hf and Th and more 
than 20% of Nb, La and U. The FLUME source region, although generating 
highly enriched magmas contains only about 10 to 20% of the Earth's K, 
Rb, Sr, Nd, Sm, Yb, Th, and U* It Is Che primary repository for Ba and 
La and has small total abundances of Y, Zr, Nb, and Hf (less than 7%). 
The crust and the two upper mantle reservoirs account for about 75% of 
the heat production from K, Th, and U# Therefore, only 25% of Che 
Earth's heat flow comes from the lower mantle# About 40% Is due to 
MORBS# The continental crust and PLUME each contribute about 18%. 

This distribution of heat sources affects the style of convection 
in the two reservoirs# Assummlng PLUME to be the shallower reservoir it 
is primarily heated from below and will therefore be characterized by 
narrow ascending plumes# MORBS Is primarily heated from within and will 
be characterized by broader ascending regions and narrow descending Jets 
or "slabs". 
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The wlnerslogy of the two source regions 

Assuming that the two source regions are i tinted » an apparent 
partition coefficient can be determined by forming the ratio of the 
concentrations in MORB to those in continental tholeiltes (CFB)> This 
is shown in Figure 2 along with mineral/melt partition coefficients for 
garnet I clinopyroxene, orthopyroxene and olivine. Tholeiltes are used 
in this comparison since they presumably are products of large degrees 
of partial melting. They will, therefore, have the trace element 
pattern of their source and a relatively uniform enrichment. The ratio 
of the concentrations in the two reservoirs can be explained If the MORE 
source Is composed mainly of garnet and clinopyroxene and PLUME has been 
enriched with fluids that were in equilibrium with this layer and 
therefore depleted in such elements as Zr, Nb, Y, Yb, and Hf. 

MORB is enriched with respect to continental tholeiltes in Y and Zr 
and only slightly depleted in Hf and Yb. Both MORB and continental 
tholeiltes represent rather large degrees of partial melting of their 
respective source region. A garnet peridotite would be expected to lose 
most of its garnet and clinopyroxene under these circumstances. The 
above results indicate that primary garnet must be a minor mineral In 
PLUME but a major component of the >!ORB source. MORB is depleted In the 
most incompatible elements, such as K, Rb, Ba, and U but not as much as 
would be expected If it were composed entirely of residual crystals. 
This suggests that MORBS evolved from a melt, presumably the result of 
an early dlfferentatlon event, and consists of a mixture of melt with 
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txcess garnet which has settled Into It during crystallization* i>e. it 
is an orthocuoulate* The upper part o£ the ciantle would consist of 
lighter* cuaulate and residual* crystals and late stage fluids which 
inpart an inverse garnet signature to this region. 

This sequence of events* based on trace elements, is precisely what 
O'Hara et al. (l4) proposed on the basis of petrological and major 
element considerations. They concluded that the parent magma for ocean 
island tholeiites had experienced a previous history of eclogite 
extraction. The remaining magma evolved to ocean island tholelite by 
olivine fractionation. On the basis of major element (20) and trace 
element chemistry, I suggest that the eclogite cumulates constitute the 
source region for mid-ocean ridge basalts. The parent magma from which 
both MORBS and PLUME evolved by crystal fractionation would be picritlc, 
the result of extensive partial melting of a primitive garnet peridotlte 
mantle. Similar considerations have led to the concept of a magma ocean 
on the moon (41). Because of the higher pressures in the Earth's mantle 
a deep eclogite cumulate is the analog of the floating plagioclase 
cumulate that forms the lunar highland crust. 

In order to see If this holds up quantitatively for the trace 
elements we investigate the following model. The primitive mantle is 
split into two reservoirs by 15% partial melting and melt extraction U, 
the surface. The melt fraction reservoir is enriched by about a factor 
of 6 in the incompatible elements but, because of the low partition 
coefficients* has nearly primitive ratios of Rb/Sr and Sm/Nd. 
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Crystslllzation of th« melt yields an eclogite cumulate layer which Is 
modelled as a 50:50 mix of garnet and melt, proportions appropriate for 
an orthocumulatc* Part of the final 5% melt fraction, a melt that Is In 
equilibrium with eclogite. Is removed to deplete the cumulate layer and 
enrich the overlying mantle* The mixing ratios were adjusted so that 
the former can yield MORB by 15% olivine fractionation and the latter 
can yield CFB by 20% partial melting (38). 

Results of this model are shown as solid points In Figure 2. They 
are consistent with MORB being the result of nearly complete melting of 
a source region which represents high pressure cumulates from the melt 
fraction of the primary differentiation, l.e* eclogite or garnet 
pyroxenite* 

Prior to providing depleted MORB this region lost Its late stage 
fluids or, alternatively, crystallized completely and then lost Its 
early melt fraction on a subsequent reheating cycle. PLUME was enriched 
by upward migration of this fluid. Thus, two stages of dlf ferentatlon 
and fractionation are required In order to concentrate the LIL In the 
two upper mantle reservoirs. Isotoplcally , the MORB source evolves as 
nearly primitive mantle until It transfers Its Incompatible elements to 
PLUME. It subsequently evolves as a depleted reservoir. Basalts from 
PLUME will have tlme-lntegrated depleted, enriched or "primitive" 
Isotopic ratios, depending on when the enrichment occurred. 
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Loot tong of th« two sourca rtgloni 

Thors to voriety of evidence that suggest! that the plume source 
region Is shollov* Anomalously low seismic velocities can be traced to 
a depth of about 250 km under Yellowstone but deeper velocities appear 
to be normal (17). The variable and small amounts of partial melting 
required to generate raelllitltes, nephelenltes, basanites and alkali 
basalts (6) from a common source region are consistent with diaplrs 
rising from various shallow depths. Tholeiltes are the result of larger 
degrees of melting, consistent with adiabatic ascent from deeper levels. 
The xenoliths entrained In kimberlites and alkali basalts are volatile 
and trace element rich and these are samples from the upper 200 km of 
the mantle (12). Intuitively, we expect that volatiles will migrate 
upward and be trapped by the cold upper mantle. Volatile and LIL 
enriched magmas occur not only under continents but also at Island arcs 
where the downgoing slab perturbs the upper mantle. Plume type basalts 
occur In continental rifts. Tliese rifts evolve to oceanic ridges, with 
Isolated oceanic Island hotspots, when the rifting has led to the 
formation of an ocean basin. 

There are also some suggestive geometric constraints. Tlie volume 
of the low-velocity zone (LVZ) Is adequate to provide the plume basalts 
but not the voluminous MORE (10). The preferred mode of convection In a 
region with a high temperature gradient and a rapidly varying viscosity 
takes the form of hexagonal cells with upwelllng centers (18). This 
seems to be consistent with hotspot patterns (19) and the small areal 
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extent of hotspots. In more homogeneous regions of the mantle, where 
the temperature gradient is smaller and the viscosity more uniform, 
linear rolls arc a possible mode of heat transport. Linear ridge 
systems and the uniformity of oceanic tholelltes therefore suggest a 
deeper source, one that Is below the large vertical and lateral 
variations that occur In the upper 200 km. The sink of oceanic 
lithosphere. Judging from the depth distribution of earthquakes, appears 
to be between 200 km and 670 km. The distinctive Isotopic and trace 
element signature of MORBS can be malntairied If It Is also between these 
depths. Upper mantle temperatures are closest to the melting point 
between about ISO and 250 km depth and this Is therefore the depth range 
where It is most likely for dlaplrs to originate. 

It Is difficult to estimate the depth at which partial melting 
first occurs. Seismic data from several hotspots give low velocities 
down to at least 150 km (17,21). Eruption temperatures of 1300*C and 
large degrees of partial melting require Initiation of melting below 
150-200 km (22). There Is some evidence for crystal fractionation in 
melts as deep as 280 km (23). Mantle fragments brought up by alkali 
basalts and kimberlites do not, in general, have the trace element 
pattern required for the source region of MORB. Therefore, PLUME magmas 
may originate from shallow depths but MORBS must be deeper than 200 km. 
In a gravitationally stratified mantle one would expect the deeper layer 
to be denser and therefore rich In garnet. The inability of young 
buoyant lithosphere to subduct below about 220 km suggests a density 
Increase at this depth and this may be the boundary between the two 
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source regions (20). Taking all evidence Into account It appears that 
PLUME nay be coincident with the LVZ and dlaplrs rising from the 
transition region supply the magmas that evolve to MORB. 

A terrestrial magma ocean 

The complementary LIL element patterns of the major terrestrial 
magma types Is reminiscent of lunar data which has led to the widely 
accepted concept of crystal fractionation and cumulate formation In a 
magma ocean or vast lava lakes (41). The various source regions are 
attributed to cumulate and residual fluid layers that resulted from 
crystallization of a '^300 km thick magma ocean that, In turn, was 
derived from very early melting of at least half of the moon. If a body 
as small as the moon experienced such extensive differentiation and 
fractionation, then the Earth should have as well. The amount of 
partial melting required, ~15%, to explain the enrichment of the upper 
mantle reservoirs on Earth, Is, In fact, relatively modest f Although 
the energy of accretion of the Earth Is much greater than that of the 
moon, the greater size of the Earth results in high Rayleigh number 
convection and rapid Increase of the melting curve with depth. This 
plus the latent heat buffer and the high melting temperature of olivine 
may prevent more extensive melting on a global scale. A 15% melt 
Implies a magma ocean ~400 km deep. The pressures In the Earth are 
greater than on the moon. Early ecloglte fractionation at depth 
therefore preempts the extensive plagloclase fractionation that resulted 
In the early lunar anorthosltlc crust. 
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Crystslllcation of a oagoa ocean will proceed froa the base because 
of the relative slopes of the adlabat and the llquldus. The near 
llquldus phases at depths greater than 60 kn are garnet and 
cllnopyroxene (14, 16) and they will forn an ecloglte or garnet 
pyroxenlte ctimulate layer* Although there Is a seismic discontinuity In 
the mantle near 400 km, ecloglte is denser than residual garnet 
peridotlte to depths of 670 km (20). These cumulates will therefore 
sink to this depth, displacing residual mantle upwards* The shallower 
part of the mantle will therefore consist of olivine and pyroxene 
cumulates ami residual fluid from the magma ocean and displaced residual 
lower mantle. Tlie top of the ecloglte layer could be as shallow as 
220 ka (20) or as deep as 370 km based on considerations of this paper. 
The mantle discontinuity near 400 km, which Is usually attributed to the 
olivine-spinel phase change could therefore be either a chemical 
discontinuity or could represent the completion of the garnet pyroxenlte 
(ecloglte) to garnetlte (garnet solid-solution) phase change (42)* 

Isotopic ratios of the two source regions 

The study of Sr and Nd Isotopes places important time constraints 
on the evolution of mantle reservoirs (24-26). The continental crust 
and the depleted reservoir have mean ages of 1*5 b*y. (3). Kimberlites 
and continental flood basalts, although enriched in LIL, have been 
attributed to a primitive reservoir (2, 3,, 24, 27)* This Interpretation 
Is not required by the data as we show later. Some mantle samples have 
been derived from ancient enriched reservoirs (28)* Other samples come 
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froa rtfsrvoirt which appcrtntly h«v<i been enriched only recently 
(28,29). 

Coneider a prlnitlve mantle that partially melta and aaparates into 
two reaervoiri, aa before, by upward removal of the melt* The lower 
mantle consists of residual crystals and is therefore depleted in the 
incompatible elements. The melt fractionates into a deep ecloglte 
cumulate layer and a shallow olivlne-orthopyroxene cumulate layer A 3% 
melt fraction is transferred from the deeper to the shallower layer at 
various times (40). The Nd and Sr isotopic ratios for the two 
reservoirs ere shown in Figure 3. The central horizontal scale gives 
the ages of the depletion and enrichment eveots. Tl^e theoretical 
Cjjd-tsr correlation lines agree with the data and Indicate that the 
mid**Qceen ridge basalt source region was depleted at times between 1.3 
and 2.3 X 10^ years ago. A redistribution of LIL makes it possible to 
satisfy the mantle isotopic data even if the primary differentiation 
occurred early In the history of the Earth. Thus, the trace element, 
petrological and thermal constraints on the evolution of the mantle and 
the various reservoirs are not contradicted by the isotopic data. The 
type of model Investigated here la similar to previous ideas of mantle 
metasomatism and trace element redistribution (30). 

Implications for mantle evolution 

Isotopic studies indicate that the two major mantle reservoirs have 
been Isolated for more than I or 2 b.y. (2,3). The evidence presented 
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h«rt from trmc* alomontB, and pravloualy from major alamanc and aaiamlc 
conaldaraclons ( 20 ), auggaat that the mantle la chanlcally stratified 
and that the various regions of the mantle are eomplamantary products of 
tarrastrlal dlffarantlatlon. The present emphasis has been on magmas 
from these reservoirs but similar conclusions result from the study of 
solid fragments from the mantle (14, 20). O’Hara et al. (14) pointed 
out that the restricted number of phases found In blmlnerallc eclogltes 
from kimberlites Indicate that they are either crystal accumulates or 
crystalline residue developed In contact with a liquid. This is exactly 
the situation inferred from the trace element patterns for the source 
region of abyssal tholelltes. Eclogite and garnet perldotlte xenollths 
from kimberlites may represent samples from MORBS and PLUME, respec- 
tively. They give satisfactory average mantle compositions for the 
major oxides when combined In the portions Indicated by the trace 
elements (20). 

Although whole mantle convection, in the conventional sense. Is 
precluded In a cliemlcally stratified mantle, transfer of material Into 
and out of the various reservoirs Is possible because of the large 
volume change associated with partial melting and phase changes, e.g., 
basalt-ecloglte. The oceanic part of the plate tectonic cycle may be 
suimnarized as follows. Partial melting in the eclogite layer allows 
dlapirs to rise to the base of the oceanic lithosphere. Nearly complete 
melting occurs during adiabatic ascent. This Is possible because of the 
proximity, in temperature, of the liquldus and the solidus In eclogite. 
Perldotlte dlapirs can only partially malt because of their high 
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llqulcll. Th« ttfultting mIc frac&lonatat in naac surface nagma chaisbers 
to form thoXclitlc melts which ara light enough to rise to the surface* 
The pyroxanita residua forms the lower oceanic lithosphere* As the 
lithosphere cools and thickens the lower part transforms to garnet 
pyroxanita or eclogite which is denser than the underlying mantle* The 
oceanic lithosphere becomes gravitationally unstable and it returns to 
its source region which lies between about 220 or deeper and 670 km 
depth* Partial melting in the shallow enriched perldotite layer 
generates continental and ocean island basalts and a harzburgite 
residue, both of which are lighter chan their source region* These 
PLUME products remain in the cruSt ard upper mantle* Some volatiles are 
returned to PLUME by subducted sediments and hydrothermally altered 
oceanic crust. The high temperature gradient In the thermal boundary 
layer at the PLUME-MORBS interface brings temperatures there close to 
the melting point (20). When mature oceanic lithosphere passes over a 
hotspot it is quite likely chat the lower pyrpxenite part will 
contribute xenoliths to the alkalic magmas and, by remelting, contribute 
to the trace element and isotopic signatures of ocean island basalts* 
The evolution of the mantle is shown schematically In Figure 4. 

It has generally been assumed chat basalts of all kinds represent 
partial melts of peridotltes. Mantle compositions based on this premise 
have much lower abundances of SIO 2 , AI 2 O 3 , and CaO chan models based on 
cosmochemlcal considerations. The evidence used to construct the 
petrological models is obtained from the upper 2C0 km of the mantle and 
the models are therefore strictly only valid for the uppermost mantle* 
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A thick ccloglcc layer sarvas to Ineraaaa tha abundancaa of tha abova 
componants and can raconcila tha coaDochanical and patrologlcal 
Intarpratatlona. Tha poaaabillty of an acloglta layar in tha uppar 
■antla waa algo diacuaaad by Praia (43), 

Tha idea that acloglta nay ba tha lourca for oceanic basalts ‘it an 
old one (31) but has not bean in favor in recant years. The main 
objection is that limited partial malting of aclogita does not generate 
a tholaiita. Extensive or complete malting is required and this has 
bean thought to be unlikely, 

Hoifivar» aclogita has an extremely small melting interval (about 
60*G) compared to paridotita (16, 32), Since diapirs cannot rise out of 
an acloglta layer into a less dense peridotite layer until they are 
already extensively molten, it requires only a small additional 
temperature rise, relative to the liquldus, to complete the melting. 
This can be accomplished in a rising diapir. It is still unknown why 
melting Initiates in the first place. One possibility is insulation by 
the thick continental lithosphere. 

Hotspots and ridges in the Atlantic and Indian oceans were beneath 
continental lithosphere prior to 200 m.y. ago. A larg«i number of 
hotspots are presently under the relatively stationary continent of 
Africa, This suggests that continental insulation, which prevents 
mantle heat from being efficiently removed to the surface, may be the 
cause of the partial melting that ultimately results In ridge and 
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hottpot volcanlpnii Tha thicki ISO km* contlnantal llehoaphara (20) 
pravanta dlapirlc uprise and excanslva smiting nay ba possible until Che 
eontinant rifts and rapid spraadlng and haat ranoval can taka placa# It 
la not so obvious chat this axplanation holds for tha rldgas and 
hotspots in tha Pacific. Tha thick crust undar tha ocaanic plataaus In 
tha central Pacific or a Chick ocaanic llchospharai dua to a period of 
slow spraadlng, could also serve to insulate the taantla« raise 
tanparacurcs and causa extensive In-sltu raaltlng prior to dlapirlc 
ascent. Parts of tha Pacific rim continents such as iVlaska, Mexico, 
Central America and S.E Asia may also have been located In the the 
central Pacific prior to 200 m.y. ago. Once Initiated, tha rise of 
upper mantle dlapirs should be rapid (33), possibly rapid enough to 
avoid crystal fractionation en route. 

Conclusions 


Garnet and possibly cllnopyroxene control the complementary trace 
element patterns of MORE and plume basalts. 1 suggest that the source 
region for MORB was formed from an ecloglte cumulate layer resulting 
from crystallization of the plcrltlc melt fraction of the original 
differentiation of the Earth. The plume source Is enriched in those 
elements chat would be concentrated in a melt that was removed from a 
garnet rich region. This fluid, Che result of an early partial melt or 
a lace stage fluid from crystallizing ecloglte in the MORB source, 
infiltrated the shallow mantle plume source at various times. This 
gives plume basalts an Inverse garnet trace element signature relative 
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to MORB* The MORB end plume source regions together make up about 
25-30X of the mantle. This requires that most, if not all, of the 
mantle has been processsed to obtain the observed enrichment. There is 
no need to invoke a large primitive reservoir in the mantle or deep, 
lower mantle plume sources. On the contrary, the plume source region 
appears to be shallow and may be coincident with the LVZ. The upper 
mantle transition region, 220-670 km, represents about 21% of the mantle 
and the LVZ is about one-third as large. This is about the ratio of the 
sizes of the two upper mantle reservoirs which is required to satisfy 
the trace element data. The MORB source is probably an eclogite 
cumulate while the plume source is probably garnet perldotite. 
Crystallization of a magma ocean would give this kind of upper mantle 
stratigraphy. Nd and Sr isotopic data indicates that the depletion of 
that part of the MORB source region being sampled today occurred over 
the Interval 1.5 to 2.5 x 10^ years ago. The fluid involved in the 
depletion/enrichment events, and in the formation of the continental 
crust, appears to be a partial melt or late stage residual fluid from an 
eclogite cumulate layer. The lower mantle is residual perldotite and, 
although it may be, convectlng, it no longer communicates with the 


shallower reservoirs 
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Figure CapClone 

Figure !• NorTiiellzed trace clcncnt concentrations In the continental 
crust (dots), continental basalts and mid-ocean ridge basalts 
(2, 6, 8, 3A, 35)* All concentrations are nortMllzed to 
terrestrial values of Canapathy and Anders (37), recalculated 
to mantle equivalents. T\ie lower curves for continental 
tholelites and MORB are normalized to mantle concentrations 
after removal of the continental crust, as are the curves for 
alkali basalts and kimberlites. 'Ihe complementary nature of 
MORB and the other magma types Is evident. The continental 
crust is not the only "enriched" reservoir. 

Figure 2. Crystal/llquid partition coefficients for garnet, 
clinopyroxrne, orthopyroxene and olivine (6). Heavy line is 
ratio of concentrations in primitive mid-ocean ridge basalts 
(36) and continental tholelites (2, 6). This Indicates that 
MORB results from melting of a garnet and clinopyroxene rich 
source region and that the continental tholeiite (CFB) source 
region has experienced ecloglte extraction. The points are 
for a model (38) in which the proto-MORB source region is an 
ecloglte orthocumulate from the initial melt fraction of 
primitive mantle differentiation. The continental and ocean 
Island basalt, or plume, source region Is depleted residue 
enriched by a late stage melt from the proto-MORB cumulate 
layer. 
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Flgurt 3. Prcssnt valuts (39) for an eclogice cuoulate 

(•olid line) which has been depleted at various times by 
removal of a melt fraction (40) « This melt fraction 
enriches the complimentary reservoir (dashed line)* Data 
from references 2, 24, 2S, 29, and 30. 

Figure 4. A model for the evolution of the mantle. Primitive mantle 
(1) Is partially molten either during accretion or by 
subsequent whole mantle convection which brings the entire 
mantle across the solidus at shallow depths. LIL elements 
are concentrated In the melt. The deep magma ocean (2) 
fractionates into a thin plagloclase**rich surface layer and 
deeper ollvlne-*rlch and garnet-rich cumulate layers (3). 
Late stage melts In the ecJoglte cumulate are removed (4) 
to form the continental crust (c.c.), enrich the perldotlte 
layer and deplete MORBS, the source region of oceanic crust 
(o.c.) and lower oceanic lithosphere. Partial melting of 
PLUME (5) generates continental flood basalts (CFB), ocean 
Island basalts (OIB) and other enriched magmas, leaving a 
depleted residue (harzburglte) that stays In the upper 
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IV. Evidence from Geophysics 

The earth is the best studied planet and It Is reasonable to test 
cosmological models against the structure and composition of the earth 
as deduced from seismology or even to proceed in the opposite direction. 

The earth appears to b <3 a differentiated body with the crust, mantle, 
core and inner core being the main subdivisions. Presently forming oceanic 
crust can be understood as a differentiated, or partial melt fraction, 
of the upper mantle and petrological models of the upper mantle agree 
well with the selsmologlcal data (1, 2). There is also no great difficulty 
in understanding most other crustal rocks in terms of a similar origin 
involving either primitive melt's or melting of reworked crustal material. 

The older massive anorthosite complexes may, however, involve processes 
or materials that were unique to the earliest history of the earth. 

The seismic data for the upper mantle (above 400 km) Is well explained 
by a mixture of olivine, pyroxene and garnet In the proportions suggested 
from petrological studies. The seismic discontinuity at 400 km can be 
explained In terms of the ollvlne-P-splnel phase change In a predominately 
olivine upper mantle, and there Is no need to Invoke any changes In com- 
position down to at least 600 km. The general Increase In seismic velocity 
and density Is well explained by self -compression and by fairly well 
understood phase changes In the ollvlne-pyroxene-garnet system. The olivine 
successively transforms to the 0-splnel and then to the spinel structure. 
Involving a total density Increase of about 10%. 

V The pyroxene and garnet components Interact to form a complex garnet 

I 

{ 
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■olid solution so that at 650 km the mantle Is composed of about equal 
parts of spinel and garnet. 

The 650 km discontinuity In the mantle represents a large Increase 
In both density and seismic velocity. It cannot be attributed to any 
known phase change in the olivine component (12) and must therefore Involve 
the pyroxenc~gamet component If It Is to be Interpreted In terms of a 
phase change alone. There is evidence » however, that the lower mantle 
does not have the same composition as the upper mantle and that the 650 
km discontinuity represents a chemical, as well as a phase, boundary (1-7). 

The absence of earthquakes below 650 km can be Interpreted In terms of a 
density barrier, due to a change In chemistry, although other Interpretations 
arc possible. The seismic data for the lower mantle is consistent with 
SIO 2 and FeO enrichment, compared to the upper mantle, and with a composition 
close to ordinary chondrites ( 1 ) . 

The lowermost mantle has unusual seismic properties and appears to 
differ in composition from the bulk of the lower mantle. It can be ex- 
plained In terms of AI 2 O 3 and CaO enrichment (7) which Implies enrichment 
in U and Th, providing a mechanism for driving the core dynamo by differential 
heating from above. 

The outer core la molten and has the properties of Iron with about 10% 
of a lighter element or compound which is probably S or 0. The Inner core 
Is solid and, within the uncertainties of present data, can have the same 
composition as the outer core or be pure Fe or Fe-Ni. 

The relative amount of iron in the earth is significantly greater 
than chondritic abundances. It is also greater than inferred for the Moon, 
Mars and possibly Venus and less than Mercury. In terms of the condensation 
hypothesis Iron is a refractory. The earth Is then presumably also enriched 
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In ocher refractories such as Ca-Al silicates and oxides and U and Th. 

The lowermost mantle may still be enriched in these refractories. 

The bulk of the lower mantle appears to be closer to pyroxene chan to 
olivine stoichiometry and closer to chondrltlc than upper mantle coiitpoaltlont 
at least in terms of the Mg, Fe, and Si ratios (1). This Is an Important 
boundary condition which is not generally considered in cosmochemlcal or 
petrological modeling. The bulk compositions of the terrestrial planets 
probably have higher Fe/Mg and Si/Mg ratios than the present terrestrial 
upper mantle which is probably a differentiate of the lower mantle. Further- 
more, from thermal and energy considerations this differentiation is likely 
to be more efficient on Che earth. The source region for basalts on 
smaller bodies is likely also to be enriched in Fe and Si relative to the 
source region of terrestrial basalts. In addition, the source regions of 
Che refractory planets. Mercury and the Moon, are probably enriched in Ca 
and A1 relative to both the earth and chondrites. 

The origin of the earth's core is difficult to understand in the frame- 
work of Che homogeneous accretion hypothesis. However, Che central concen- 
tration of U and Th and possibly ^^Al makes it easy to melt the core and 
maintain the dynamo and the associated refractories can explain the lower- 
most mantle (8, 9, 13). 
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Plan«e«ry Modelling 

Although the earth can be used as a guide in constructing models of 
the terrestrial planets the above differences should be kept in mind, i.e.» 

X) the earth is enriched in iron and possibly other refractories compared 
to other inner solar system objects with the exception of Mercury and 2) 
the upper mantle of the earth is not representative of the average com- 
position of the mantle. The phase changes in a purely pyroxene planet 
are quite different from those in a mainly olivine planet. 

As a first cut* however* the density depth profiles* and hence the 
mean density and moment of inertia of the terrestrial planets can be con- 
structed from the terrestrial 'density-pressure curve. Thermal models can 
be used to correct for differences In temperature. This affects not only the 
density versus depth curves but also the locations of the temperature 
sensitive phase boundaries. This can have a large effect on the Inferred 
mantle composition. 

This approach was used by Anderson (10) on Mars and it was concluded 
that the Martian mantle was enriched in FeO compared to the terrestrial 
upper mantle and that a generally satisfactory model could be constructed 
with a chrondritic starting material. The moment of inertia Implies that 
Mars has a small core. The data can be satisfied by an Incompletely 
differentiated carbonaceous chrondrite (minus most of the water) or a 
mixture of chrondrltes. 

When the iron content of a planetary mantle is Increased the location 
of the equilibrium phase boundary between olivine and the post-olivine 
phases shifts to lower pressures or shallower in the planet. The effects 
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of composition and temparature on the location of phase boundaries are 
■ore important than refinements in the equations. of state. In small 
planets, such as the Moon and Mara , temperature is as important as pressure 
in controlling the density but temperature, of course, is much more uncertain 
than pressure. 

Pressure vs. depth is almost an invariant in planetary models that 
have the same radius, mass and moment of Inertia. It can be computed 
from the assumption of hydrostatlclty. Even if the planet is grossly out 
of hydrostatic equilibrium the effective pressure, which controls the’ 
density, can be estimated quite accurately. The actual density at a given 
depth is another matter. There is a trade-off between temperature and 
composition. VIhen we have only two parameters, such as mean density and 
moment of inertia, we can only determine two parameters and these can be 
for example, two of the following three: density of the mantle, density 

of the core and radius of the core. If information about the crust is 
available, from gravity or seismology, this can be stripped off. 

For small planets, where pressure is not great, a very simple equation 
of state can be used, such as the Murnighan or Birch-Murnlghan or, as 
discussed above, the ’’terrestrial equation-of-state." The zero-pressure 
density in these equations is a variable, to be determined by the mass 
and moment of inertia, or a parameter set by assumptions about composition. 

A simple equation of state and one which is adequate for small planets 
is the Murnaghan equation, 

P - ^ l(p/Po)" - 11 


> ■ 
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where Is the zero-pressure bulk nodulus» p Is the density » Pq is the 
zero pressure density end n is the pressure derivative of the bulk nodulus, 
n ■ (dK/dP)o» i^ich is generally from 3 to 4. To a first approxlaetion 
the effect of teapereture can be includid simply by writing, 

P - <k(T)/n) t<p/po(T)*^ - 13 

A higher order equation is the Blrch-Mumighan equation which is a 
■saber of a class of equations which can be written in general form, 

p - 3Ko (.-nJ-lKp/Po)^^**"*^’- 

where m is generally 1 or 2 and n is 2 or 4 with n < n. Again, to first 
order the effect of temperature can be included in and Pq. Equivalent 
expressions can be written for the seisnlc velocities (11). 

The pressure is, 

P(z) • / p(z) g(z) dz 

where, 

and z and r are depth and radius and g Is the gravity. 

Much more sophisticated equations have been applied to the earth 
but they depend on a detailed knowledge of the seismic velocities which 
are unknown on the other planets. 

For the moon the effect of pressure is small and a direct comparison 
of seismic velocities can be made with only a small correction for pressure. 
Uncertainties in the temperature and temperature coefficients of the physical 
parameters are the main sources of uncertainty in constructing and inter- 
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ThTwml History Modtlllng 


In order to calculate the thermal evolution of a planet the following 
■uat be known: 

1) The temperature and velocity of the Infalllng material, 

2) The rate of Inf all and the energy redistribution upon Impact 
and embedment, 

3) The ratio of retalnment and loss of Impacting material as a function 
of size, volatility, velocity and time, 

A) The thermal and elastic properties of the material as a function 
of temperature and pressure, 

5) The surface boundary .condition and opacity of the environment 
as a function of time, 

6) The radioactive abundances, half-lives and absolute time of 
initiation of accretion. Short-lived radiogens may be Important If accretion 
Is early and rapid, 

7) The mode of cooling, l.e. conduction, radiation, convection, 
and 8) The latent heats of all phase changes Involved in Impact and 
compression processes. 

Clearly, approximations and assumptions must be made. The end member 
assumptions are a) that all kinetic energy Is retained In the body and 
Used up as specific heat and latent heat and b) the body grows In radiative 
equilibrium with the surroundings. The efficiency with which a body 
retains energy brought In depends on how deeply the particles are embedded, 
how much material Is vaporized or ejected and how rapidly they are covered 
up. Thus, crater mechanics and ballistics are Involved and It Is Important 
to know the size and Impact velocity of the accreting material. Thus, a 
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amount InConoatlon la requlrad bafora th& ^calculation can ba 
■tartad and sona of thaaa dapand on coaaolo|y» i.a, ate tha partlclaa 

coming from and «hat haa baan thalr prevloua hlatory? 

Ttia quaatlon of homogeneoua veraus Inhomoganaoua atcratlon is Important 
at the earllaat atagas. Whan vaporization and aacapa are Important tha 
planet retalna only tha refractor lea and "homoganaoua" accretion actually 
leads to an Inhomoganeous planet. Whether a planat starta off aa rafraetory 
depends not only on tha temperature In tha nebula and the composition of 
the planatasifflals but on where In the solar system the particles are 
originating. Advocates of pure homogeneous accretion will have to demonstrate 
that Incoming planetealmala either escape vaporization or, If rot, that the 
volatiles can be retained by ths planet. These conditions can be expected 
to Vary with time. 

The present internal temperatureo of a planet are under a little 
bit better control since the planet strives to evolve to an equilibrium 
configuration. A hot planet will simply approach equilibrium faster than 
a cold planet. The present internal temperatures si'.i.ply represent a 
balance between heat production, conduction and convection, the latter two 
being a function of the initial temperatures. The boundary conditions are 
the surface temperature and if the gravity field is known, an estimate 
of the thickness of the cold outer boundary layer. The age and depth of 
origin of magmas is a potential constraint. High thermal gradients can 
be expected in the outer layers and a near adiabatic gradient at deeper 
layers. The temperature of the bulk of the interior is therefore set by 
the thickness and conductivity of the outer shell. 
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